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We investigate the high-order behavior of perturbative matching conditions 
in effective field theories. These series are typically badly divergent, and 
are not Borel summable due to infrared and ultraviolet renormalons which 

i-^H ■ introduce ambiguities in defining the sum of the series. We argue that, when 

treated consistently, there is no physical significance to these ambiguities. 

K^ . Although nonperturbative matrix elements and matching conditions are in 

H I general ambiguous, the ambiguity in any physical observable is always higher 

order in 1/M than the theory has been defined. We discuss the implications 
for the recently noticed infrared renormalon in the pole mass of a heavy quark. 
We show that a ratio of form factors in exclusive A;, decays (which is related to 
the pole mass) is free from renormalon ambiguities regardless of the mass used 
as the expansion parameter of HQET. The renormalon ambiguities also cancel 
in inclusive heavy hadron decays. Finally, we demonstrate the cancellation 
of renormalons in a four- Fermi effective theory obtained by integrating out a 
heavy colored scalar. 
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I. INTRODUCTION 

In physical problems involving several distinct scales, it is often convenient to describe the 
physics using an effective field theory. Typically, one is interested in physics at energies much 
less than the mass of a heavy particle, in which case the physics is most easily described by 
an effective Lagrangian in which virtual heavy particle exchange is accounted for through a 
series of non-renormalizable operators. The coefficients of these operators are perturbatively 
calculable as a power series in as{M), where M is the mass of the heavy particle which is 
integrated out. However, since the resulting perturbation series is only asymptotic, and 
since it is well known that perturbative QCD is not Borel-summable, the perturbatively 
calculated coefficient functions are at some level ambiguous. One source of this ambiguity 
is infrared renormalons ||l]-§l, which arise in QCD from graphs of the form shown in Fig. 1. 
They contribute to the factorial growth of the coefficients of the perturbation series, and 
introduce uncertainties in the sum of the perturbation series proportional to powers of 
Aqcd/M. 

There has been much recent discussion in the literature on the effects of infrared renor- 
malons in the context of one particular effective field theory, the heavy quark effective theory 
(HQET) |§-0] . In particular, the use of the pole mass as an expansion parameter in HQET 
has been criticized, because it has been shown to suffer from an ambiguity which prevents 
its definition to an accuracy better than C(Aqcd) [1,11 , and the authors of Ref. p advocate 
formulating HQET in terms of a short- distance mass m(/i) which does not suffer from a 
renormalon ambiguity at 0{Aqcd)- The results of Ref. |^,^ clarify the formal status of the 
pole mass (and of nonperturbative matrix elements in an effective theory in general). 

Given this ambiguity, it is important to show that physical predictions of HQET, which 
are usually expressed using the pole mass as an expansion parameter, are unambiguous. 
In this paper, we investigate the effects of renormalons on matching conditions in effective 
field theories. We argue that, while perturbatively calculated coefficient functions suffer 
from renormalons, any ambiguity in a physical observable is always higher order in 1/M 



than the theory has been defined and is consequently irrelevant. Therefore, as long as 
one works consistently, it does not matter that unobservable parameters such as the heavy 
quark mass or the matrix elements of higher dimension operators are not unambiguously 
defined; relations between physical quantities are unambiguous. We also argue that, while 
formulating HQET in terms of some short-distance mass m{fi) is certainly possible, use of 
an expansion parameter other than the pole mass (or some mass which differs from rripoie by 
0{Aqcd)) is both inconvenient and uneccessary. 

The paper is arranged as follows: In Section 2, we review the physics of infrared renor- 
malons, and discuss the basic method of calculation that will be used in this paper. In 
Section 3 we discuss HQET and the renormalon ambiguity in the quark pole mass. We 
show explicitly that there are no renormalons in the ratio of form factors for A5 -^ Ac 
semileptonic decay from which the perturbatively defined pole mass may be extracted, and 
we comment on the renormalon cancellation for inclusive semileptonic B meson decay. In 
Section 4, we discuss renormalon ambiguities in a four-Fermi effective theory. We show that 
this effective theory has new renormalon ambiguities not present in the full theory which 
cancel corresponding ambiguities in matching conditions, and also show that the cancella- 
tion of renormalon ambiguities is not specific to HQET, but occurs in other effective field 
theories. This section can be omitted by readers only interested in HQET. The conclusions 
are presented in Section 5. 

II. RENORMALONS 

QCD perturbation theory is used to express some quantity / as a power series in a^, 

00 
fias) = /(O) + J2U «r^ (2.1) 

n=0 

Typically, this perturbation series for / is only asymptotically convergent. The convergence 
can be improved by defining the Borel transform of /, 

00 r 

B[m) = f{0)6{t) + Y.^t-, (2.2) 

71=0 "• 



which is more convergent than the original expansion eq. ( p.l|) . The original expression 
f{as) can be recovered from the Borel transform B[f]{t) by the inverse Borel transform 

/(a.)= / dte-*/"=i?[/](t). (2.3) 

Jo 

If the integral in eq. ( |2.3| ) exists, the perturbation series /(os) is Borel-summable, and is 
unambiguously defined. However, if there are singularities in B[f]{t) along the path of 
integration, the function / is ambiguous. The inverse Borel transform must be defined 
by deforming the contour of integration away from the singularity, and the inverse Borel 
transform in general depends on the deformation used. 

One source of singularities in B[f] in QCD is infrared renormalons ^-^. These arise 
from graphs of the form in Fig. 1. Physically, these graphs correspond to the running of a^, 
and infrared renormalons are ambiguities in perturbation theory arising from the fact that 
the gluon coupling gets strong for soft gluons in the one- loop diagram in Fig. 1. The infrared 
renormalons produce a factorial growth in the coefficients /„, which gives rise to poles in 
the Borel transform B[f]. The renormalon ambiguities have a power law dependence on 
the momentum transfer Q^. For example, a simple pole at t = to in B[f] introduces an 
ambiguity in / depending on whether the integration contour is deformed to pass above or 
below the renormalon pole. The difference between the two choices is proportional to 

5/~ / e-*«/"=(«) 5[/](t) (2.4) 

^[ Q ) 

where bo = —{11 — 2r7,j/3)/47r is the leading term in the QCD /3-function 

^"^2 = boa's + O {a!) , (2.5) 

which governs the high energy behavior of the QCD coupling constant 



and the contour C encloses to- It is useful to write the Borel transform B[f]{t) in terms of 
the variable u = — fco^-O The form of the renormalon singularity in eq. ( p.4| ) then implies 
that a renormalon at uq produces an ambiguity in / that is of order (Aqcd/Q) ""■ This 
ambiguity is canceled by a corresponding ambiguity in non-perturbative effects such as in 
the matrix elements of higher dimension operators. The sum of the perturbation series plus 
non-perturbative corrections is expected to be well defined. 

The Calculational Method 

Clearly, one cannot sum the entire QCD perturbation series to determine the renormalon 
singularities. Typically, one sums bubble chains of the form given in Fig. 1 0,§]. Beneke § 
considered a limiting case of QCD in which the bubble chain sum is the leading contribution 
to the renormalon. Take QCD with Nf flavors in the limit Nf —>■ oo with a = NfOg held 
fixed. Feynman diagrams are computed to leading order in a^, but to all orders in a. Terms 
in the bubble sum of Fig. 1 with any number of bubbles are equally important in this limit, 
since each additional fermion loop contributes a factor asNf, which is not small. QCD is 
not an asymptotically free theory in the Nf -^ oo limit, so the procedure used by Beneke 
is to write the Borel transform as a function oi u = —h^t but still study renormalons for 
positive u. The singularities in u are taken to be the renormalons for asymptotically free 
QCD. This procedure is a formal way of doing the bubble chain sum, while neglecting other 
diagrams. 

The Borel transform of the sum of Feynman graphs containing a single bubble chain can 
be readily obtained by performing the Borel transform before doing the loop integral [§,§. 
The bubble chain sum is 

G(a., fc) = f; (^ - g,\ {boasNjr {H-k'/f^') + C)" (2.7) 



*This is the definition used in l6|, and is the negative of the definition used in 



where k is the inomentum flowing through the gauge boson propagator, C is a constant that 
depends on the particular subtractions scheme, and Bq = 1/6tt is the contribution of a single 
fermion to the /5-function. In the MS scheme, C = —5/3. The Borel transform of eq. ( p.7| ) 
with respect to agNf is {u = —hot) 



1 1 (KK \ 



2 

2+n [^fJ.'^i' ~ ^ 9nu 



Mln(-Fe^//i') (2.8) 



1 h^y 1 i... _P 
asNf \e^ j (-P) 



The Borel transformed loop graphs can be computed by using the propagator in eq. ( p.8|) 
instead of the usual gauge boson propagator 

[k^k, - k^g,,) -^. (2.9) 



III. RENORMALONS IN THE HEAVY QUARK EFFECTIVE THEORY 

A. Matching Conditions 

An effective fleld theory Lagrangian is an expansion in an operator series in inverse 
powers of some mass scale M. By construction, the effective fleld theory has the same 
infrared physics as the full theory. However, because the ultraviolet physics (above the scale 
at which the theories are matched) differs dramatically in the two theories, the coefficients 
of operators in the effective theory must be modified at each order in as{M) to ensure that 
physical predictions are the same in the two theories. 

Since the two theories coincide in the infrared, these matching conditions depend in gen- 
eral only on ultraviolet physics and should be independent of any infrared physics, including 
infrared renormalons. However, in a mass-independent renormalization scheme such as di- 
mensional regularization with MS, such a sharp separation of scales cannot be achieved. 
It is easy to see how infrared renormalons creep into matching conditions. Consider the 
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familiar case of integrating out a W boson and matching onto a four-Fermi interaction (we 
will discuss a variant of this example in detail in Section 4). The matching conditions at one 
loop involve subtracting one-loop scattering amplitudes calculated in the full and effective 
theories, as indicated in Fig. 2. For simplicity, neglect all external momenta and particle 
masses, and consider the region of loop integration when the gluon is soft. When k = 0, 
the two theories are identical and the graphs in the two theories are identical. This is the 
well-known statement that infrared divergences cancel in matching conditions. However, 
for finite (but small) k, the two theories differ at 0{k'^/M^^) when one retains only the 
lowest dimension operators in the effective theory. Therefore, the matching conditions are 
sensitive to soft gluons at this order, and it is not surprising that (as we shall show) the 
resulting perturbation series is not Borel-summable and has renormalon ambiguities starting 
at OiAlcD/M^). 

However, this ambiguity is completely spurious, and does not mean that the effective 
field theory is not well defined. Since the theory has only been defined to a fixed order, an 
ambiguity at higher order in l/M^r is irrelevant. The renormalon ambiguity corresponded 
to the fact that the two theories differed in the infrared at (9(A;^/M^). When operators sup- 
pressed by an additional power of 1/M^ in the effective theory are consistently taken into 
account, the two theories will coincide in the infrared up to (9(A;'^/M^), and any ambiguity 
is then pushed up to 0{AQu^/M^r)■ Consistently including l/M^r suppressed operators 
pushes the renormalon to 0{AQfj^/Myy), and so on. In general, a renormalon at -u = -Uq in 
the coefficient function of a dimension D operator is canceled exactly by a corresponding 
ambiguity in matrix elements of operators of dimension D — 4 + 2uo, so that physical quan- 
tities are unambiguous. This cancellation is a generic feature of all effective field theories, 
and also occurs in HQET. 

B. HQET and the Quark Mass 

The HQET Lagrangian has an expansion in inverse powers of the heavy quark mass. 



-^0 = hy (^D ■ v)hy — 5m h^h^, (3.1) 

£i = . . . . 

Here ^ught is the QCD Lagrangian for the hght quarks and gluons, h^ is the heavy quark field, 
and Ck are terms in the effective Lagrangian for the heavy quark that are of order I/ttIq. 
There are two mass parameters for the heavy quark in eq. ( p.lD , the expansion parameter 
of HQET rriQ, and the residual mass term 6m. The two parameters are not independent; 
one can make the redefinition mo —>■ wt-q + Am, 6m, — > 6m, — Am. A particularly convenient 
choice is to adjust m-o so that the residual mass term 6m, vanishes. Most HQET calculations 
have been done with this choice of mg, but it is easy to show that the same results are 



obtained with a different choice of mo ||T0[ . The HQET mass when 6m, = is often referred 
to in the literature as the pole mass, and we will follow this practice here. 

Like all effective Lagrangians, the HQET Lagrangian is non-renormalizable, so a specific 
regularization prescription must be included as part of the definition of the effective theory. 
An effective field theory is used to compute physical quantities in a systematic expansion 
in a small parameter, and the effective Lagrangian is expanded in this small parameter. 
The expansion parameter of the HQET is Aqcd/'^o- One can then use "power counting" 
to determine what terms in the effective theory are relevant to a given order in the 1/mo 
expansion. For example, to second order in 1/mo, one needs to study processes to first 
order in C2, and to second order in Ci. It is important that the renormalization procedure 
preserve the the power counting for the effective field theory to make sense. In order to 
preserve power counting, we must choose a mass-independent subtraction scheme; in our 
case, we choose to use dimensional regularization and MS. A mass-dependent subtraction 
scheme (such as a momentum space cutoff) mixes different orders in the 1/mo expansion. 
Thus to compute a quantity to first order in 1/mo, one would have to retain the effective 
Lagrangian to all orders in l/m,Q, which is not particularly useful. Hence, HQET is defined 
using a mass-independent subtraction scheme, and nonperturbative matrix elements must 
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be interpreted in this scheme. 

It has recently been shown [^,0 that there is a renormalon in the relation between the 
renormalized mass at short distances (such as the MS mass m) and the pole mass of the 
heavy quark at u = 1/2, which produces an ambiguity in the relation between the pole mass 
and the MS mass of order Aqcd- This implies that there is an ambiguity in the residual 
mass term 6m of order Aqcd due to renormalon effects |^,||. 

The quark mass in HQET and the MS mass at short-distances are parameters in the 
Lagrangian that must be determined from experiment. Any scheme can be used to com- 
pute physical processes, though one scheme might be more advantageous for a particular 
computation. The MS mass at short distances is useful in computing high-energy processes. 
However, there is no advantage to using the "short distance" mass (such as the running 
MS mass) in HQET, as advocated by [Q. In fact, from the point of view of HQET, this is 
extremely inconvenient. The effective Lagrangian eq. ( p.l|) is an expansion in inverse powers 
of mQ. Power counting in l/rriQ in the effective theory is only valid if 6m is of order one (or 
smaller) in tuq, i.e. only if 6m, remains finite in the infinite mass limit m,Q -^ oo. When m,Q 
is chosen to be the MS mass the residual mass term 6m is of order mo (up to logarithms), 
which spoils the l/mg power counting of HQET, mixes the a^ and 1/mo expansions, and 
breaks the heavy flavor symmetry. For example, using rrio to be the MS mass at /i = ttiq, 
one finds at one loop that 

4 

6m = --a^mo. (3.2) 

671 

In 6 — i> c decays, including this residual mass term in the heavy c-quark Lagrangian, causes 
l/m,c operators such as hc{—'i^)Th[,/mc to produce effects that are of the same order in 1/mc 
as lower dimension operators of the form hcThi,. While physical quantities calculated in this 
way must be the same as those calculated using the pole mass, it unnecessarily complicates 
calculations to use a definition for ttiq that leaves a residual mass term that is not finite 
in the m,o — » oo limit. Better choices of the expansion parameter of HQET are the heavy 
meson mass (with 6m, of order Aqcd), and the pole mass (with 6m, = 0). 



The MS mass at short-distances can be determined (in principle) from experiment with- 
out any renormalon ambiguities proportional to Aqcb/itt-q (i-e. at m = 1/2). As an example, 
consider the MS mass of the 6-quark at the GUT scale in an SU{5) unified field theory. The 
6-quark mass at the GUT scale is proportional to the 6-quark Yukawa coupling at the GUT 
scale, which in turn is equal to the r-lepton Yukawa coupling at the same scaleQ There 
are no QCD renormalons at u = 1/2 in the relation between the r-lepton mass at short- 
distances, and the pole mass of the r (neglecting QED effects). Thus one could determine 
the 6-quark mass at short distances by measuring the r-lepton mass, without any renormalon 
ambiguities at m = 1/2. 

The MS quark mass can be related to other definitions of the quark mass using QCD 
perturbation theory. The connection between the Borel transformed pole mass and a short- 
distance mass (such as the MS mass) has been worked out in ref. 0. The relation between 
the two (for the c-quark) is 



T71 



37rNf 



^J -^ ^^^--) r(3-u) -^ + ^^^(-) 



(3.3) 



where rric is the renormalized [not the pole) mass at short distances, such as the MS mass, 
/x is the renormalization scale, and the constant C and the function i?Si (w) depend on the 
renormalization scheme. Eq. ( ^.31 ) has a renormalon singularity at m = 1/2 which is the 
leading infrared renormalon in the pole mass. Writing u = 1/2 + Am, we have 

B[ml^'%u = 1/2 + An) = --^ — + . . . (3.4) 

STTi\fmc Am 

where the ellipsis denotes terms regular at Am = 0. In the next two sections, we will only 
work to (9(l/r?7,o), so poles to the right of m = 1/2, which are related to ambiguities at 
higher order in l/rriQ, are irrelevant at this stage. 



T^There are corrections to this relation from matching conditions at the GUT scale, which will 
have renormalon ambiguities proportional to powers of Aqcd/"^gut- 
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Although m^°^^ is formally ambiguous at 0{l/mc), we will argue in this paper that phys- 
ical quantities which depend on mP°^'^ are unambiguously predicted in HQET. We demon- 
strate this explicitly for a ratio of form factors in A;, semileptonic decay. We then comment 
on the cancellation of renormalon ambiguities in the expression for the inclusive semileptonic 
width of the B meson. Both results will make use of eq. (|3.4|) and its analog for the b-quark. 



C. A from Exclusive Decays 

The matrix element of the vector current for the semileptonic decay A;, -^ AcC" z/e decay 
is parameterized by the three decay form factors 

{Ac{v')\cYb\Ab{v)) = u{v') [Fi{v ■ v')Y + ^2(1; ■ v')v^ + F^iv ■ v')v'^] u{v). (3.5) 

In the limit ruh, rric — > 00, and at lowest order in a^, the form-factors F2 and F3 vanish. We 
will consider a^ and l/rric corrections, but work in the nib = 00 limit. Consider the ratio 
rp = F2/F1, which vanishes at lowest order in a^ and l/rric- The corrections to rp can be 
written in the form [^ 

rF{as,v-v) = -— - = — -— - + fr{as,v-v) 3.6 

Fi[v ■ v'j mc[l + V ■ v') 

where the function f^{as,v ■ v') is a perturbatively calculable matching condition from the 
theory above fi = rric to the effective theory below fi = rric, and the A term arises from l/rric 
suppressed operators in HQET. At one loop, |n]] 



/,(«„ vv') = -^ , \^ In {vv'+^iv v'f - 1 ) . (3.7) 

"J^ \ [v ■ v'Y — 1 



The ratio rp = F2/F1 is an experimentally measured quantity, and does not have a 
renormalon ambiguity. The standard form for rp in eq. ( ^.6| ) is obtained by using HQET 
with the pole mass as the expansion parameter. The HQET parameter A is the meson mass 
in the effective theory, i.e. it is the meson mass mp, minus the pole mass of the c-quark. The 
pole mass has the leading renormalon ambiguity P,^ at m = 1/2 given in eq. (|3.4| ), which 
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produces an ambiguity in the l/rric contribution to F^jFi given by the first term in eq. ( ^.61) . 
There must therefore also be a renormalon at u = 1/2 in the radiative correction to F2/F1 
given by the second term in eq. ( |3.6| ). It is straightforward to show, using the techniques of 
Section 2, that this is indeed the case. 

The Borel transformed series B[f^]{u,v ■ v') in the 1/A^/ expansion is easily calculated 
from the graph in Fig. 3 using the Borel transformed propagator in eq. ( p. 81) . The Borel 
transform of the Feynman diagram is 



The radiative correction to F2 (which determines /r) is obtained from the terms in eq. ( p.8| 
which are proportional to f ". Combining denominators using the identities 

1 , . fi dxil-xY^'' 



{2 + u)j 
Jo 



nS+u' 



(fc2 + 2meA;-t;')(A;2)2+" ' 'Jo [k^ + 2m^xk-v 

1 / X /""^ dX 

2(3 + m) 



{k^ + 2m^xk ■ vf^"" k ■ V Jo [k^ + 2m^xk ■ v' + 2Xk ■ v]'^^''' 

extracting the terms proportional to f '^ and performing the momentum integral, we obtain 



3TcNf{l + u) \e^ J Jo Jo [X"^ + 2XmcXV ■ v' + mlx'^] 



B|MK''-'0- J:",, :^,, (^) ."./, 'iA(dx ,„,„,;;^^;;_, ,--_^,^„.,., (3.9) 



Rescaling A —>■ xrricX and performing the x integral gives 



3TTNf\mle^'J T{3-u) Jo [A2 + 2At; ■ t;' + 1]'+"' 

(3.10) 

This expression has a pole at m = 1/2. Expanding in Am = u — 1/2 gives 



^ L rdx - 

37rmce^/2 Am io [A^ + 2At; • f ' + 1]' 
2/i 1 1 



B[fr]{u = 1/2 + Au,v v') = ,__^,, — I dX ^_ _ —73J2+- 



(3.11) 



STcrrice'^/'^ Am 1 + m ■ m' 
where the ellipsis denotes terms that are regular at n = 1/2. 

The Borel singularity in eq. ( pJ.llD cancels the singularity in the first term of eq. ( ^.61 ) at 



all values of v ■ v', so that the form factor ratio rF{as,v ■ v') = F2{v ■ v')/Fi{v ■ v') has no 
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renorinalon ambiguities. Therefore the standard HQET computation of the l/mc correction 
to F2/ Fi using the pole mass and the standard definition of A gives an unambiguous physical 
prediction for the ratio of form factors. 

D. Inclusive Decays 

A similar situation occurs for inclusive B decays, which have been the subject of much 



recent interest ||12|-|1^ . The inclusive B -^ Xgcu (where q = u or c) decay rate is related to 



the imaginary part of the forward scattering amplitude, 

T{B -^ Xgeiy) ~ Im (5|T(J'^^ J'-')\B), (3.12) 

where J^ = C7'^(l — 75)6. In this case the expression for the total rate as an expansion in 
powers of l/mq is not the result of matching onto an effective theory, but instead is the 
result of performing an operator product expansion on the time ordered product of the two 
currents in Eq. ( ^.121 ). The final expression is 



TiB^X.eu) = ^l|^|L!!^[/o(a,) {B{v)\ hh \B{v)) 

+ —fiias) {Biv)\ hiiD-v)h \Biv)) + Oiml/ml,AlaD/ml)], (3.13) 
nil. 



where /o = 1 and f\ = 1 to lowest order in a^. Eq. (|3.13|) is true with an arbitrary residual 



mass term in the HQET Lagrangian, and we have not yet applied the equations of motion 
to the operator hbi{D ■ v) hi,. The total decay rate F is an observable, and does not have a 
renormalon ambiguity. It was shown in Refs. 0,^ that the total decay rate is unambiguous. 
It is important to note that this result does not require the use of a "short-distance" mass 
in Eq. ( p.l3| ). One is free to choose some definite (but arbitrary) prescription for integrating 
around the pole at u = 1/2 in Eq. ( |3.3| ). The mass rrih in the leading term of Eq. ( |3.13| ) 
is then well-defined, but there is an ambiguity at (9(Aqcd) in the residual mass term Sm 
arising from the renormalon at m = 1/2 in the pole mass. By the equations of motion, 

[{v-B) h = 5mh + 0{l/mb), (3.14) 
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so the ambiguity at (9(Aqcd) in the matrix element {B{v)\ hbiD ■ vhb \B{v)) is the same as 
that in the pole mass, eq. ( p.4| ). From ||^, the Borel transformed series B[fQ{u)] for the a^ 
corrections to the leading term has a pole at u = 1/2, 

i;|/„|(.= l/2 + A,,) = ^^^^£^^ + ... (3,15) 



Comparing eq. ( |3.4] ) with eq. ( |3.15| ) and eq. ( |3.13| ), we see explicitly that the ambiguity in 



the matrix element oi iD ■ v cancels that in /q to give an unambiguous prediction for the 
total width F. 

Thus in this O.P.E. the cancelation of renormalon ambiguities occurs in the same manner 
as in the construction of HQET: the ambiguity in the matrix element of a higher dimension 
operator cancels that in the perturbation series for the coefficient of the leading operator. 

IV. FOUR-FERMI THEORY 

As we argued in Section 3, the cancellation of renormalon ambiguities between matrix 
elements and matching conditions is a general feature of an effective field theory. In this 
section, we illustrate this in a more familiar effective field theory, four- Fermi theory in which 
a heavy colored scalar is integrated out. (We choose the scalar theory as our example because 
it has Feynman graphs which are slightly easier to compute than in the four-Fermi theory 
of weak interactions.) The theory is QCD with Nf light fiavors in the Nf —>■ oo limit, with 
agNf fixed. Two of the Nf fiavors (called d and s) couple via a color triplet scalar of mass 
M ^ Aqcd to color singlet particles (called e and r) according to the interaction Lagrangian 

C = A (sr + Je) + h.c. . (4.1) 

We will study the theory in the d ^ s sector. 

The effective Lagrangian obtained from eq. (|4.1| ) after integrating out the heavy scalar 
has the form 

A^ A^ A^ 

^«// = ^ ^^(^) ^^ ^^ + j^^ ^^(/^) ^^/^^^ ^ ^^''^ ~ ^ ^^(^) ST D'^[ed] + ... (4.2) 
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where the elhpsis indicates additional operators of dimension 6 and higher. The scale /x is 
the scale at which the effective Lagrangian is computed from the full theory, and is usually 
chosen to be /i = M . The effective Lagrangian at zeroth order in the strong interactions is 



computed by equating the dr -^ se scattering amplitude obtained from eq. ( |4.2| ) with that 
obtained by expanding the scalar exchange graph in the full theory (Fig. 4) in a power series 
in 1/M^, to give cs = I, ct = 0, and ds = 1. 

We now consider matching onto the effective theory at high orders in as, and concentrate 
on the coefficient of the operator Ot = sa^^r ea^'^d. This is a particularly convenient 
operator because its coefficient ct is zero at tree level, and at one loop in the full theory 
it only receives a contribution from the graph in Fig. 5. The other graphs do not have the 
correct 7-structure to contribute to ct, and contribute only to the scalar amplitude, ct is 
obtained by equating the tensor scattering amplitudes in the full and effective theories, and 
its Borel transform is computed using the techniques discussed in Section 2. 

The Borel transform of the tensor scattering amplitude A^ (at zero momentum transfer) 
in the full theory is obtained by evaluating the Feynman graph in Fig. 5, 



B 



.full 



' ^ Usa^^Ur Uea^fsUa / -—7 — ^^^TT^ „.,,, ,,2+« (4-3) 



using the gluon propagator Eq. ( ^.8] ). The Ui are Dirac spinors for the external fermion 

lines, and the quarks have been given a common mass m to regulate the infrared behavior 

of the diagram. Performing the k integral gives 

■ \ 2 / 2 \ " 



B 



full 



/\rp 



l27iN V^j ^'>''"^^r^-''^^-^dT{u-l)T{2-u) (4.4) 

" (M^)'~" - (m^)'" [(1 -u)M^ + urn''] 
_ (M2 - m'^f 

Despite the appearance of V{u — 1) in eq. ( [4.3| ), the amplitude is finite at u = and 1 as the 

term in the square brackets vanishes at these points, but it has singularities at m = 2, 3, . . .. 

However, we stress that renormalon singularities in Eq. ( |4.4| ) are of no interest, since we 

are not going to attempt to calculate scattering amplitudes in perturbation theory. Any 

attempt to do so will of course face serious infrared problems. We are only interested in 
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using perturbation theory to calculate the coefficient functions in Eq. ( [4.2|) , which requires 
us to subtract the corresponding amplitude calculated in the effective theory. 
To match onto the effective theory, we expand (|4.4|) in a power series in 1/M : 



• \ 2 / 2 \ " 

B[A'f] = j^ I ^ U^a'^-'Ur Ue<y,uUd T{u - I) 1(2 - u) (4.5) 



/ 



X 



(-T" ,_,,(:^,... 



M2 ' ' M2 

where we have only retained terms up to order 1/M^, since ct/M"^ is the coefficient of a 
dimension six operator. It is perhaps useful at this point to relate this to the standard 
perturbation series for A^^\ 

^^"" = \2'kN^M^ ^^^''^^ ^^^^'^^'^ ^ «n(«siV/)" + 0{—). (4.6) 

Taking the appropriate derivatives of ( |4.5| ) gives for the ffist few terms of the series 

Oo = 0, 

ai = -l-21og(m/M), 

as = -6o [C + 2(1 + C) log(m/M) - 2 log(/i/M) 

+21og^(m/M) - 41og(m/M) log(/^/M)] , 
as = ... (4.7) 

The tensor scattering amplitude in the effective theory is computed from the loop cor- 
rection to the lowest order operator c^, 



B 






1671 X'^ f fi^Y - ^^ - ^^ f d^k k'^k^ik'^k'' -eg"") 



^Nf \e^J -^-^^--^----^-'^ J (2vr)4 (p _ ^2)2 ^2 (_fc2)2+«' 

where we have used the tree-level value C5 = 1 in evaluating the graph. This result is the 
same as that obtained by setting the scalar propagator in eq. ( |4.3| ) to l/M^, since the cs 
term in the effective Lagrangian reproduces this piece of the four- Fermi vertex. Evaluating 
the k integral gives 

B[Af] = ^^^ l^^j U^a^-^U^ U^a,Ma r(^) r(2 - ..) i^. (4.9) 
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This reproduces the (m^) term in eq. (|4.5|) , including the entire u dependence. Comparing 
eq. (^75|) with eq. (|4.9|), we obtain 



^1=-(''>1 = I2^(?) nu-l)r(2-u)LJ_, (4.10) 



Note that any dependence on the quark mass m has dropped out of ( |4.1CI| ), so that the 



matching condition is independent of the infrared regulator. Therefore, in terms of the 
original perturbation series 

\2 CO 

CT(/i) = -^^ Y^CniasN^r, (4.11) 

all large logarithms of m/M have dropped out of the matching conditionsj^ 

Co = 0, 

Ci = C"-21ogWM), 

C2 = -bo [C" + 2Clog(/i/M) - 2 log' (/i/M)] , 

C3 = ... (4.12) 

(where C and C" are scheme- dependent renormalization constants). However, despite the 
fact that the individual terms q in the expansion of ct are now well-defined, the expression 
( [4.10| ) has poles at u = 0, 1, 2, .... The pole at m = is removed by renormalization 0, but the 
renormalons at m = 1, 2, ... correspond to ambiguities of order (Aqcd/^)'" in the coefficient 
function ct{^i)- Note that these are different from the singularities in the expression ( [4.41 ). 
The singularity at m = 1 in ( |4.10| ) is not present in ( [4.4| ), while the coefficients of the 



singularities at m = 2, ... also differ. For the singularity at m = n we find 



•fThe logarithms of m/M are reproduced in the effective theory by scaling the operators from 
/i = M to low energies. 
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where the elhpses denote terms regular at u = n. Note that the singular piece of i?[A^^'] at 
u = n has no term proportional to (yU^/M^)". 

We now consider using the effective Lagrangian in Eq. ( |4.2| ) to calculate the cross-section 
for the spin-flip scattering process rd -^ es. This is of course not a physical process, since 
the d and s quarks are not physical asymptotic states, and the perturbatively calculated 
rate A^^^ exhibits serious infrared problems, as is clear from Eq. ( [4 .41) . However, we may use 
this process to demonstrate that physical predictions in the effective theory are well-defined, 
by demonstrating that the renormalon ambiguities cancel between the coefficient function 
cxifJ') and the (perturbatively computed) matrix elements of higher dimension operators. 
Since the graphs which cancel the ambiguity also occur for the matrix elements of physical 
hadrons, this cancellation will also take place for physical amplitudes such as r — > K^e. 

The order 1/M^ contribution to rd -^ se process from the operator Ot is 

' ^ ^ T{u-1) r(2 - u) ^-—^ U^a^-'Ur U^a^Md. (4.14) 



Writing u = 1 + An, one finds that the singularity at m = 1 is 

^' ^41 774 f^^^"'^- ^ecr,.t/. -^ (4.15) 



which is of order 1/M^. The renormalons at w = 2, . . . produce singularities of order 
1/M^, . . .. Since we have only computed the effective Lagrangian to order 1/M^, we can 
ignore the renormalon singularities at u > 2, and the only singularity that is relevant to 
the order we are working is the one at u = 1. This singularity is canceled by a singularity 
in the Borel-transformed matrix element of the ds operator, which is also of order l/M''. 
The matrix element of the sr D^ [ ed ] operator between quark states is evaluated using the 
graphs of Fig. 6, where only the first graph contributes to the spin-flip scattering amplitude 
l^TiX^ i^^Y- ,, - ,, ( d^k k'^k^ikf'k'' -eg'"')k'^ 



Evaluating the k integral gives 

'^' ^^yU,a^'^UrU^a,M,nu-l)Ti3-u) "^^ ^fj \ (4.17) 



UnNfye^'J ' ^r^e^>.u^a^,^ ^; ^ v^ ^, ^4 



Expanding around u = 1 gives the singular term 

BlA^^h = ^ 1^1 Usa^^'Ur UeOauUd 4t ^- (4.18) 



This is precisely the negative of eq. ( [4.15| ), so that the singularity cancels in the total 



amplitude, which is the sum of the two terms eq. ([4.14| , [4.17|) . 



Clearly, at m = 2, 3, ... similar cancellations will take place with 1/M^, 1/M^, ... operators. 
This is simply because the singularities in B[ct{^i)\ are not found in the scattering amplitudes 
in the full theory, Eq. ( |4.13|) . Since the full and effective theories are, by construction, 
identical up to the order to which the effective theory has been defined, the singularities 
must cancel between matching conditions and the matrix elements of higher dimension 
operators in the effective theory, as they do at u = 1. 

V. CONCLUSIONS 

Perturbatively calculated matching conditions in an effective field theory suffer from 
renormalon ambiguities. However, we have argued that any ambiguity in a physical quan- 
tity is always higher order in 1/M than the effective theory has been defined and is there- 
fore of no consequence. In practice, one calculates matching conditions to a given number 
of loops, and from physical measurements then determines the value of nonperturbative 
matrix elements. The cancellation of renormalon ambiguities in physical observables then 
simply means that, although the values obtained for the nonperturbative matrix elements 
will depend sensitively on the number of loops at which the theories are matched, relations 
between physical quantities will not. If the unphysical parameters are extracted from ob- 
servables at a given order in a^, then they can be used to predict other observables to the 



same order in a^, as was done for example for the extraction of |Vbc| in |T^,[T^. Inclusive and 
exclusive semileptonic decays of hadrons containing a heavy quark are free of renormalon 
ambiguities, regardless of the mass parameter of the l/irtQ expansion. In addition, we have 
demonstrated renormalon cancellation in an effective field theory other than HQET — the 
four-Fermi theory. 
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Figure Captions 

1 The bubble chain diagrams which are the leading contribution to the renormalon as 
the number of light colored fermions is large, Nf -^ oo. 

2 The one-loop contribution to the matching of a higher dimension operator. The coef- 
ficient of the operator in the effective Lagrangian is cq + ciag + C2Q;^ + ... • This series 
will have IR renormalons due to the incomplete cancellation of the soft gluons in the 
two graphs. 

3 The loop graph with the Borel transformed gluon propagator contributing to the v'^ 
form factor in A^ -^ AcCu. 

4 The tree- level exchange of a heavy colored scalar. It only contributes to cs{n) and 
dsifJ') and not to Ct{^)- 

5 The leading contribution to ct{h) arises at one loop in QCD. The blob on the gluon 
propagator denotes the bubble sum of light quark loops. The external quarks have 
been given a small mass m which regulates the IR behavior of the graph. 

6 The contribution to the scattering amplitude from the {iDY operator which is 
0{1/M^). Only the first diagram contributes to the spin-fiip amplitude at order 1/Nf. 
The ambiguity at -u = 1 in this matrix element cancels that in Ct{^i)- 
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Figure 1. 
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Figure 2. 




Figure 3. 
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Figure 4. 




Figure 5. 
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Figure 6. 
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